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Within a truncated determinant algorithm, two alternatives are discussed for including systematically the 
remaining ultraviolet modes. Evidence is presented that these modes are accurately described by an effective 
action involving only small Wilson loops. 
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1. Introduction 

Because QCD in four dimensions is renormal- 
izable, not super-renormalizable, the fluctuations 
of the fermion determinant are significant at all 
physical scales. Fortunately the short distance 
behaviour of QCD is very well understood. In 
particular, we know that for sufficiently high mo- 
mentum scales this physics should be accurately 
described by an improved gauge action. 

In the truncated determinant algorithm pi, the 
fermion determinant is separated into two pieces 



lndet/f = [TrlniT 



low X 



[Trlni?]^ 



gh A 



(1) 



where the lowest n cut eigenvalues are directly cal- 
culated and included in the Monte Carlo updat- 
ing procedure. The higher eigenvalues can be in- 
cluded in the Monte Carlo by some approxima- 
tion that matches onto the low eigenvalue results 
without gaps or double counting, is controlled and 
becomes exact in the continuum limit. 

Two numerical methods suggest themselves for 
calculating the high eigenvalues: (l)The multi- 
boson approach of Liischer||. (2)Using a small 
number of gauge loops to model the determinant 
as proposed by Sexton and Weingarten and 
Irving and Sexton B. 



2. Matching onto the Multiboson Method 

One method to compute the high eigenvalues 
which is guaranteed to succeed is the multiboson 
approach of Liischer||. Define 

P eff (U) = [det(D + m)pexp(-S s ([/)) (2) 
and 

H = 75 (D + m)/[c m (8 + m)] (c m > 1) (3) 

where c m is chosen so that the eigenvalues of H 
are in the interval (—1, 1). Liischer chooses a se- 
quence of polynomials P n {s) of even degree n such 
that 

lim Pjs) = lis for all < s < 1 (4) 

n — >oc 

then for rif = 2 

detiJ 2 = lim [detP^H 2 )]- 1 (5) 

n — >oo 

Choose polynomials such that complex roots 
Z\ . . . z„ come in complex conjugate pairs (non 
real) so that \fz = fi + iv. Then 



dettf 2 = lim TT det[(iJ - [i k f + v 2 k \ 



(6) 



fe=i 

Hence we can write 
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P ef f(U)= lim i- / D<f>D<f>Uxp-(S g + S b ) (7) 
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where the bosonic action is given by 

a 

s b = ]T ]T \{h - fi k )Mx)\ 2 + » 2 \<t>kM 2 - (8) 



fc=l x 



Liischer used Chebyshev polynomials to estimate 
how many boson fields (rib) are required to repre- 
sent the original action to a fixed accuracy in the 
range (e < s < 1) The error is given by|Q]: 



ii?( S )i<2(i-^r +i . 

1 + \/e 



(9) 



Therefore, the convergence is exponential with 
rate as rib —> oo. 

A practical problem with this multiboson 
method is that it requires an increasingly large 
number of boson fields as the quark mass be- 
comes lighter. As m q — > 0, we must take e — * 0, 
but to obtain a fixed level of accuracy we must 
hold 2^/erib fixed and hence rib increases without 
bound. 

However the multiboson method matches 
nicely onto the calculation of low eigenvalues. 
This was first suggested by Alexandrou et.al. Js| . 
In the truncated determinant method, the cutoff 
e for the multiboson method is set by the highest 
eigenvalue of H 2 which is explicitly included in 
the low end calculations. Hence it does not ex- 
plode as the quark mass goes to zero. The combi- 
nation of methods remain accurate for all quark 
masses. For example, for j3 = 5.9 on a 12 3 x24 lat- 
tice with direct inclusion of the lowest 100 eigen- 
values, the associated cutoff for the multiboson 
simulation of the high eigenvalues is y/e w 0.035 
independent of the light quark mass. 

Furthermore, the error associated with the in- 
accurate behaviour of the polynomial fit in the 
range < s < e can be corrected as low eigenval- 
ues are computed for every configuration update. 
We obtain a reweighting term, 



AS b = J2\n(\ 2 P(\ 2 )) 



(10) 



which can be included to eliminate errors in the 
region < s < e. 
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Figure 1. Comparison of (Trlni/)hi g hA (dashed 
line) and best fit (up to 6 links) effective gauge 
action (open circles) for 75 gauge configurations. 



3. Matching onto Small Loops 

Using the multiboson method for the high end 
of the determinant satisfies all our requirements 
and completes the algorithm. However it may 
be possible to reduce the total required computa- 
tions even further using a more physical approach 
to the high eigenvalues. Consider how many of 
the high eigenvalues actually have physical infor- 
mation and are not just lattice artifacts. For ex- 
ample, for a 12 3 x24 lattice with = 5.9 and K = 
.1587 there are 497,664 total eigenvalues of the 
Wilson-Dirac operator; while for a high energy 
cutoff of 1 GeV we have approximately 1500 eigen- 
values (0.3%). For a fixed volume V and quark 
mass m q a decreasing fraction of the eigenvalues 
are below a fixed physical scale as f3 — * oo. There- 
fore, most of the range of large s fit in Liischer's 
multiboson method is physically unimportant. 

This suggests a more physically motivated 
method for dealing with the high eigenvalue part 
of the fcrmion determinant in which one approx- 
imates the ultraviolet contribution to the quark 
determinant with an effective gauge action: 



[Trlntfk 



igh A 



»=0 



otiLi 



(11) 



where each Li is a set of gauge links which form 
a closed path. The natural expansion is in the 
number of links. For zero links Lq is just a con- 
stant, for four links we have a plaquette, and six 
links give the three terms found in considerations 
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of improved gauge actions |6j . 

This idea was studied in detail by Irving and 
Sexton{|. These studies were done on a 6 4 lat- 
tice at /? = 5.7 with Hybrid Monte-Carlo full QCD 
simulations (with a heavy sea quark). Their re- 
sults were rather discouraging. It was hard to get 
a good approximation to the determinant with a 
closed set of loops and they needed large loops to 
even approach a reasonable fit [3] . 

There are however two important differences 
between their study and our situation. First, they 
simulated the whole determinant, while here we 
only need to approximate the eigenvalues above 
some cutoff. Hence we would expect the small 
loops to dominate at least for sufficiently high 
cutoff. Second, they used an approximate pro- 
cedure to estimate stochastically the logarithm 
of the determinant needed, while we are exactly 
computing all eigenvalues for this study. It turns 
out that these differences are critical, as using 
approximately the same lattices (and with even 
lighter quarks) we find an excellent approxima- 
tion to the high end with only small loops. 

We generated a set of 75 configurations on a 6 4 
lattice at f3 = 5.7 and k = .1685. We included 
the lowest 30 eigenvalues (which corresponds to 
a physical cutoff of approximately ~350 Mev) in 
the Monte Carlo accept /reject step in the gener- 
ation of these independent configurations. 

Considering only the high eigenvalues, an ex- 
cellent fit to the fluctuations is obtained includ- 
ing four and six link closed loops. The variance of 
the fit is 0.265. The comparison between the fluc- 
tuations in the exact (St) and approximate (S a ) 
actions for the high eigenvalue piece is shown in 
Fig [l]. As expected, if only the plaquette term 
is included the variance is larger (2.25) and we 
must move the low eigenvalue cutoff to N = 50 
(w 700 MeV) to reduce the variance below one. 
The results for various cutoffs and terms included 
are shown in Table 1. 

The linear combination (.46,-. 55, .04, .70, .03) 
for (plaquette, rectangle, chair, polygon, wilsonlinc) 
gives the best fit to the high eigenvalues of the 
quark determinant (with n cu t = 30). The con- 
figuration to configuration variations of the indi- 
vidual 6-link terms are highly correlated. This 
is to be expected since these three terms are not 



n\ 


A 


4 links 


6 links 


cut 


(MeV) 




(with WL) 








4.98 


1.074 (0.835) 


±15 


340 


2.25 


0.2652 (0.233) 


±50 


700 


0.940 


0.0564 (0.0491) 


±250 


1,210 


0.0733 


0.0695 (0.0641) 


±1250 


2,220 


0.138 


0.0198 (0.0180) 



Table 1 
Variance 
values of 
of small e 



(< (S a - S t ) 2 > 1/2 ) of fit to high eigen- 
the quark determinant by various sets 
;auge loops (WL denotes a Wilson line) . 



independent. 

The coefficients of the effective action should 
be independent of the physical volume with other 
physical parameters held constant. We are re- 
peating this study for an (8 4 ) lattice. The prelim- 
inary study on 41 configurations gives the same fit 
parameters within the statistical accuracy. With 
n cut = 120 (A w 350 MeV) the variance is 0.52. 

Although more study is required, this second 
method looks very attractive for dealing with the 
high end of the fermion determinant in full QCD 
with light dynamical quarks. Simulations would 
be performed by including the predetermined ef- 
fective gauge action S a in the gauge updates and 
computing the infrared part of the determinant 
as in the truncated determinant simulations. 
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